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We investigate the mass of ttie parity partner N*{l/2~) of the nucleon 7V(l/2+), 
in lattice QCD using a new lattice discretization scheme for fermions, domain 
wall fermions (DWF). DWF possess exact chiral symmetry and flavor symmetry, 
both of which are required for this problem, even at finite lattice spacing. Our 
calculation reproduces the large mass splitting between those two states, in good 
agreement with experiment. We also present preliminary results for the mass of 
the positive-parity excited state Af'(l/2+). 



1 Introduction 

It is important to understand the hadron mass spectrum from the first 
principles of quantum chromodynamics (QCD), the fundamental theory of the 
strong interactions. As is well known, the only method for such a first-principle 
calculation is numerical lattice QCD. Recent lattice calculations of the light- 
hadron mass spectrum in the cuienched approximation agree well with experi- 
mental values within about 5%u. However, this success is restricted to ground 
states and does not apply to excited states. The excited-state mass spectrum 
is big challenge in lattice QCD. 

In this study, we focus on a negative-parity nucleon excited-state A^*(1535), 
which is theoretically identified as the parity partner of the nucleon A^(939). 
We are especially interested in the mass splitting between N and N* from the 
viewpoint of parity partners. As is well known, the mass splitting between 
parity partners would be absent if chiral symmetry were preserved El. In other 
words spontaneous chiral-symmetry breaking is responsible for the absence of 
such parity doubling in the actual hadron spectrum. In this sense, regardless 
of a model or a theory, it is important to handle the chiral symmetry and 
its spontaneous breaking for reproducing precisely the mass splitting between 
parity-partner hadrons. 

Indeed, both non-relativistic quark models □ and bag models □, which ex- 
plicitly break chiral symmetxy, fail to reproduce the large mass splitting be- 
tween iV(939) and -/V*(1535)El. The non-relativistic quark models are based on 
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a harmonic oscillator description of the orbital motion of constituent quarks. 
The plausible value of its oscillator quantum should be a couple of hundred 
MeV tojeproduce the observed charge radius and magnetic moment of the 
nucleonH. Since this model regards N* as a state with one quantum excitation 
in orbital motion, it indicates that N* lies just a few hundred MeV above the 



e wrong ordering of 
because A^' should 



ground state. Even worse, we face a serious problem of t 
iV*(1535) and the positive-parity excited state JV^' (1440) 
be assigned two oscillator quanta in this modelH. Furthermore this wrong or 
dering cannot be improved by the conventional one-gluon-exckange potential 
model of the residual interaction between constituent quarks In addition, 
this difficulty led GJozman and Riska to propose another candidate for the 
residual interajetionla. It is easy to see that the bag models face essentially the 
same problema. 

It is clearly an interesting question whether lattice QCD can reproduce 
this large N-N* mass splitting. However, conventional lattice fermion schemes 
had some difficulty in this challenge. The Nielsen- Ninomiya no-go theorem Q 
dictates that either chiral symmetry or flavor symmetry or both are supposed 
to be violated at finite lattice spacing, while both of them are essential in this 
subject. The Wilson fermions explicitly break chiral symmetry at finite lattice 
spacing and hence are quite inadequate for the current problem. Although 
Kogut-Susskind (KS) fermions have a remnant U{1) chiral symmetry at finite 
lattice spacing, they are still not capable of the N* mass calculation. The main 
reason is that KS fermion has only discreta.flavor symmetries which belong to 
a subgroup of the S'J7(4) flavor symmetry Q. There are only three irreducible 
representations, 8, 8' and 16 for KS baryon operators due to this incom- 
plete flavor symmetry. Two appropriate representations 8 andJ-G, to which 
iV*(1535) belongs, involve also A(1405), A(1520) and iV(1520)i. The study 
of N* spectrum using KS fermions always faces this inevitable contamination 
from lower mass states. 

Several ago, Kaplan! advocated a new type of lattice fermion scheme, 
and ShamirE3'M reformulated it for lattice QCD simulations. They are called 
domain wall fermions (DWF), which utilize a fictitious fifth dimension. An 
important feature of DWF is that chiral symmetry is almost exactly preserved 
even at the finite lattice spacing. This is achieved because the symmetry 
violation is suppressed exponentially in terms of the number of lattice sites Ng 
in the extra dimension EiJ. In other words, Ng gives us a way to control the 
violation. Recent lattice calculations with DWF have shown that good chiral 
properties are obtained for moderate like 16 if the lattice spacing is small 
enough, like 0.2 fmtil. In addition, the flavor symmetry is also well preserved 
in this fermion discretization. Thus we are led to an attempt to use DWF for 
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lattice QCD calculations of the N* mass spectrumi 



2 Baryon operators 

The mass ms of the low- lying baryon B is extracted from the two-point 
correlation function composed of the baryon interpolating operator Ob- 

Gos W = X](0|T{Os(f , t)dB{0, 0)}|0} , (1) 

which behaves like exp(— m^i) for large t since-Go^ (i) is projected out at 
zero spatial momentum through the sum over We focus on the nuclcon 

channel, the spin- half iso-doublet baryons. There are two possible interpolat- 
ing operators for the corresponding quantum number = 1/2+ even if we 
restrict them to contain no derivatives and to belong to the (i , 0) © (0, |) chiral 
multiplet under SU{2)l ® SU{2)r[A: 

= eabc{ulCj5db)uc , (2) 

B2 = £abc{ulCdb)-f5Uc , (3) 

where abc, ud, C and 75 have usual meanings as color, flavor, charge conjuga- 
tion and Dirac indices. The superscript "-f" denotes the positive parity. 

The operator is preferred for use in lattice QCD to extract the signal 
of the nucleon ground-state. On the other hand, the operator B2 is con- 
ventionally discarded since B^ is expected to couple weakly to the nucleon 
ground-state due to the vanishing non-relativistic limit 113. In our calculation, 
the nucleon interpolating operator is assigned to B^ in the conventional way. 
We try to calculate the excited-state mass spectrum from _B^. 

Multiplying the left hand side of the previous positive parily operators by 
75, we obtain the interpolating operators with negative parity Ej: 

Bi = J5B1 = eabc{ulCj5db)j5Uc , (4) 
B2 = I5B2 = £abc{ulCdb)Uc ■ (5) 

The point to notice is the relation between the nucleon two-point functions 
with opposite parities 

GB+it) = -l5GB-it)l5- (6) 

This means that the two-ppipt correlation function can couple to both positive 
and negative parity states EJ: 

GB+{t)^ (l+7*)GoW + (l-7t)Go(-t) , (7) 
Gb- (t) ^ -(1 - 7*)Go W - (1 + 7t)Go(~i) , (8) 
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where G'o(t) = 6'(t)A+e"^^+* + 6'(-t)A_e+*^-* f\. M± denote masses of the 
opposite parity state. Since Go(t) possp^es the contribution from both positive 
and negative parity lowest-lying statestJ, which propagate in the opposite time 
direction respectively, we are threatened by contamination from the backward 
propagating opposite parity state in the case of the baryon spectrum with 
some boundary condition. An appropriate boundary condition is required in 
the time direction to reduce the contamination of the opposite parity state. 

At the end of this section, let us briefly review the knowledge that Wjia,- 
broken chiral symmetry imposes parity doubling in the hadron spectrum EltJ. 
For the sake of simplicity, we consider a particular transformation of the SU{2) 
chiral symmetry as [Q5,u] = +i"f5U and [Q5,c?] = —175^. Then, one can 
easily find that the two-point correlator composed of each Bi and B2 should 
transform as 

[Qs, B+^{x)B+M = * {75, B+^{x)B+^{0)} (9) 

in the chiral limit; [Qi,,H] — 0. Suppose that the vacuum possesses chiral 
symmetry: Q5IO) = 0. Eq.(^ shows that the two-point correlation function of 
the nucleon and 75 have to anti-commute, 

{75,GB+(t)} = 0. (10) 

Immediately, Eq.(H) gives 

GBAt) = GB-{t) , (11) 

which means that parity doubling arises in the nucleon channel due to chiral 
symmetryEJ. Of course, the chiral symmetry is spontaneously broken, QsjO) 7^ 
in the QCD vacuum so that such a parity doubling never occurs in the actual 
spectrum □. In this sense, the spontaneous breaking of chiral symmetry is 
responsible for the absence of parity doubling. ^ 

3 Numerical results 

We generate quenched QCD configurations on a IG'^ x 32 lattice with the 
standard single-plaquette Wilson action at /3 = 6/5^ = 6.0. The quark prop- 
agator is calculated by using domain wall fermions with a fifth dimension of 

tin the quenched approximation, N* can be regarded as the stable baryor»4ike the nucleon. 
-fThis argument ignores the consequence of the 't Hooft anomaly condition tZI and also spoils 
the possibility of the massless baryon because Gg+ {t) and Gg - {t) are defined at zero spatial 
momentum. 
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A's^lfi, sites and the domain-wall height M^l.8. According to our simula- 
tionsEa, this corresponds to lattice units a^^ « 1.9 GeV from amp— 0A00{8) 
in the chiral limit and spatial lattice size La « 1.7 fm. 

In this work, we provide forward-type and backward-type quark propa- 
gators with the wall source at two different locations (t = 5 and 27) and 
take the average of two measurements for the hadron spectrum in each config- 
uration. We use 205 independent gauge configurations for the lightest two 
quark masses, m ~ 0.02 and 0.03 and 105 configurations for the heavier 
ones, m — 0.04 — 0.125. Those bare quark masses correspond to mass ratios 
m-^/nip « 0.59 — 0.90. All calculations were done on the 600 Gflops QCDSP 
machine at the RIKEN BNL Research Center. 

Now let us touch upon technical details. We take a linear combination of 
two quark propagators with periodic and anti-periodic boundary condition in 
the time direction, respectively. This procedure enables us to extract the state 
with desired parity even in the baryon spectrum. 

3. 1 Parity partner of nucleon: N* 

We present mass estimates of the nucleon (N) and its parity partner (N*) 
obtained by the single mass-fit method applied to the two-point functions. 
We first calculate the effective masses to find appropriate time ranges for fit- 
ting. The effective mass plot shows plateaus on the far side of the source 
(t — isourcc=13-20) for N and on the near side of the source (t — isourcc=5- 
12) for N*. We choose to fit our baryon propagators from some minimum 
time slice Tmin to an appropriate maximum time slice Tmax — 20 for N and 
Tmax < 16 for N* . To keep fitting ranges as wide as possible, Tmin is varied 
from Tmax — 2 and selected under the condition x^/^uf < 1-0 where A^df 
denotes the fitted degree of freedom. All our fits are reasonable in the sense 
that the confidence-level is larger than 0.3 and estimates from the weighted 
average of the effective mass agree with them within errors. 

In Figure 1 we show the low-lying nucleon spectrum as a function of the 
quark mass, m. The nucleon mass is extracted from the operator. We 
omit the point at m = 0.02 for the operator i?^ since a good plateau in the 
effective mass plot is absent. N* mass estimates from and B2 operators 
agree with each other within errors in the whole quark mass range, as expected 
from their common quantum numbers E3. Note that this result disagrees with 
the N* spectrum obtained in Ref 16: we cannot find any discrepancy between 
mass spectra from B^ and B2 ■ Note also the same signal for N* is obtained 
from a mixed correlation function (0\B^ B2 + -B^-Bj" |0) in our study. 

The most remarkable feature in Figure 1 is that the N-N* mass splitting 
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Figure 1: A'' (x) and TV* (□ and O) masses versus the quark mass in lattice units (a~^ fa 
1.9GeV from amp=0. 400(8) in the chiral limit). The corresponding experimental values for 
N and N* are marked with stars. The N-N* mass splitting is clearly observed. 



is observed in the whole range and even for hght valence quark inass valuesli3. 
This mass splitting grows as the valence quark mass decreases 113. To make 
this point clear, we compare two mass ratios, one from the baryon parity 
partners niN* / and the other from pseudo-scalar and vector mesons m-^/nip 
in Figure 2. Experimental points are marked with stars, corresponding to non- 
strange (left) and strange (right) sectors. In the strange sectoLwe use S and 
E(1750) as baryon parity partners and K and K* for mesons eI. We find the 
baryon mass ratio clearly grows with decreasing meson mass ratio, toward the 
experimental valuesE3. 

Finally, we evaluate the N and N* mass in the chiral limit. We simply 
take a linear extrapolation in 5 lightest quark masses for and . We find 
TOAr=0.55(l) and m^r* =0.80(2) in lattice units for_values in the chiral limit. If 
the scale is set as a"^ ~ 1.9 GeV from the p massEj, we obtain niN ~ 1.0 GeV 
and totv ~ 1-5 GeV in a good agreement with the experimental value. Above 
errors do not include systematic errors due to finite volume, lattice spacing 
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Figure 2: Mass ratio of the negative-parity excited-state and positive-parity ground-state 
baryons versus mass ratio of the pseudoscalar meson and vector meson. All calculation is 
done for three degenerate valence quarks. 



and quenching effects. Sucli a systematic study will be addressed in future 
calculations. 



3.2 Positive parity excited-state: N' 

In principle, the mass spectrum of the "first" excited nucleon can be de- 
rived from the two-mass fits for the nucleon interpolating operator . How- 
ever, large statistics are required. Also it is difficult to control the systematic 
errors. Indeed, several attempts to extract N' mass spectrum failed to repro- 
duce experiment. Here, we take another approach to N' spectrum. 

As we mentioned, the i?^ operator vanishes in the non-relativistic limit. 
Thus we expect that weakly couples to the nucleon ground-state since the 
non-relativistic description of the nucleon was quite successful in the naive 
quark model. On the other hand, nobody has succeeded in extracting the 
signal of the nucleon by using the so-called unconventional operator B2 in 
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Figure 3: Mass estimates obtained from (x), (o) and mixed type correlator (O). 
The experimental values for N and N' are marked with stars. 



lattice This suggests that B2 has neghgible overlap with the nucleon 

ground-state and might give us a signal for an excited-state nucleon. Of course, 
this naive expectation is against the common sense that interpolating opera- 
tors with the same quantum numbers should give the same mass spectrum. 
However, _ssie find different plateaus in effective mass plots from and B2 
operatorsllj. 

We apply the single-mass fit method described above to the two-point 
function composed of B^ ■ We find that masses extracted from B^ and B^ 
arc quite different as shown in Figure 3. We conclude that we can identify 
B2 with the "first" positive-parity excited-state {N') of the nucleon for heavy 
quarks (m ~ 0.07) E£l. The main reason is that we obtain mass estimates 
consistent with B^ from the mixed correlation function {OlB^B^ + B^BflO). 
This suggests {0\B+\N) « 0. Indeed, we see \{0\B+\N) /{0\B+\N')\^ < IQ-^ 
from the two-mass fits for the correlation at m = O.IO and 0.125. This 
is plausible since the operator is expected-to couple weakly to the ground 
state of the nucleon as we mentioned earlier This feature weakens in the 
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lighter quark mass region (from around m = 0.05). Unfortunately, we have no 
data of the mixed type correlation for in = 0.02 and 0.03. To make a definite 
conclusion about the result of B^, we need further calculations for the light 
quark mass. 

Finally, we want to touch upon the reason why we see a clear signal for 
the first time in this study while previous studies failed to do so. This should be 
related to mixing induced by explicit chiral symmetry breaking at finite lattice 
spacing which is absent in our calculation but severe in those calculations. 
Although and i?^ do not mix in the continuum because of different chiral 
structure^ it is known that unwanted mixing between them arises through the 
explicit breaking of chiral symmetry by conventional lattice fermionstJ. On the 
other hand, this breaking in DWF is expected to be suppressed exponentially 
with Ns- Thus DWF can drastically reduce the unwanted mixing EJ. Indeed 
what we see here suggests such reduction is significant. As a result, we are 
able to numerically confirm the expected feature (0|i3^|A^) ~ in the region 
of heavy valence-quark mass. 

4 Conclusion 

We studied the mass spectrum of the nucleon excited-states in quenched 
lattice QCD using the domain- wall fermions (DWF) which preserves the chiral 
and flavor symmetries almost exactly. Most importantly we demonstrated that 
this method is capable of calculating the excited-state mass of N* . 

We made systematic investigation of the N* spectrum by using two distinct 
interpolating operators, and B^ ■ We obtained mutually consistent results 
for the N* mass spectrum from both of them. This is in contrast with the 
positive parity case as described below. In practice B^ is preferable to B^ in 
extracting the N* mass signal. 

We found definite mass splitting between N* and N in the whole quark 
mass range studied. Furthermore, this splitting grows with decreasing quark 
mass. This is the first time such a remarkable feature has been observed in 
any theoretical calculation. The N* mass and the N-N* mass splitting in the 
chiral limit obtained by extrapolation is consistent with the experimental value 
within about 5-10%. This is very encouraging for further investigations of N* 
physics. 

We observed no signal for the nucleon by using the unconventional nucleon 
operator which vanishes in the non-relativistic limit. Instead, we extracted 

-l-To speak properly, — B2 andpSj'" + belong to two distinct (^,0) ffi (0, ^) chiral 
multiplets under SU{2)l >^ SU{2)ii[d. 
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the mass of the first excited nucleon N' for heavy quarks. We need further 
study to make a definite conclusion about this subject for fighter quarks. 
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